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ABSTRACT 


An  18  foot  tall  jacket  type  structure  was  tested  for  three  dimen¬ 
sional  forces  and  deck  deflections  using  regular  and  irregular  waves 
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of  up  to  2.5  feet  in  height.  Theoretical  studies  for  both  regular 
and  irregular  wave  forces  are  presented  in  conjunction  with  experi¬ 
mental  results.  In  addition,  a  similar  comparison  for  deck  deflec¬ 
tion  and  frequency  response  is  presented  in  this  report. 

Results  indicate  base  shear,  lift  force,  and  overturning  moment 
calculations  using  the  range  of  hydrodynamic  coefficients  recommended 
by  the  American  Petroleum  Institute  are  in  general  agreement  with 
experimentally  determined  values  for  wave  frequencies  between  0.3 
and  0.8  Hz.  Deck  deflection  measurements  were  significantly  greater 
than  those  calculated  using  standard  structural  matrix  techniques. 

The  difference  is  concluded  to  be  the  result  of  inadequate  analytical 
description  of  the  test  structure  joints  using  a  rigid  joint  analysis. 
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NOMENCLATURE 


A  wave  amplitude 

C  wave  celerity 

Cq  drag  coefficient 

Cj^j  inertia  coefficient 
D  member  diameter 

E  modulus  of  elasticity 
F  total  force 

Fq  magnitude  of  forcing  function 
H  wave  height 

H$  significant  wave  height 
I  moment  of  inertia 

K  member  stiffness 

K‘  effective  spring  constant 
L  wave  length 

M  total  moment 

M'  effective  mass 

Sp  force  spectral  density 
SM  moment  spectral  density 
Sn  wave  spectral  density 
T  wave  period 

Tp  force  transfer  function 
Tm  moment  transfer  function 
U  deck  deflection 

anx  water  particle  acceleration  in  the  x  direction  normal  to  member 

any  water  particle  acceleration  in  the  y  direction  normal  to  member 

anz  water  particle  acceleration  in  the  z  direction  normal  to  member 

ax  water  particle  acceleration  in  x  direction 

a2  water  particle  acceleration  in  z  direction 

cx  x  angle  relationship  of  member  orientation 

Cy  y  angle  relationship  of  member  orientation 

cz  z  angle  relationship  of  member  orientation 

f  force  per  unit  length 
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fm  frequency  of  peak  energy 

fx  force  per  unit  length  in  the  x  direction 

fy  force  per  unit  length  in  the  y  direction 

f2  force  per  unit  length  in  the  z  direction 

g  gravitational  constant 

h  water  depth 

k  wave  number 

i  member  length 

u  horizontal  water  particle  velocity 

urms  root  mean  square  of  horizontal  water  particle  velocity 

w  vertical  water  particle  velocity 

un  horizontal  water  particle  velocity  normal  to  member 

wn  vertical  water  particle  velocity  normal  to  member 

x  member  phasing  distance 

z  depth  below  the  still -water  level 

a  tower  orientation  measured  clockwise  from  tank  axis 

e  phase  angle 

C  damping  ratio 

n  water  level  displacement  from  still -water  level 

e  angle  from  +x  axis  to  member  axis 

\>  resultant  water  particle  velocity 

p  density  of  water 

<j>  angle  from  +z  axis  to  member  axis 

03  angular  frequency 


SIGN  CONVENTION 
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INTRODUCTION 


For  over  30  years  the  offshore  oil  industry  has  been  employing 
fixed  offshore  structures  for  use  in  the  recovery  of  natural  gas  and 
oil  reserves.  In  addition,  military  and  defense  requirements  have 
also  required  use  of  fixed  offshore  structures  for  many  years.  For 
example,  the  offshore 'Texas  Towers' of  the  early  1950s  formed  part  of 
the  early  warning  system  used  by  the  United  States  and  the  Navy's  air 
combat  maneuvering  range  off  the  coast  of  North  Carolina  utilizes 
fixed  offshore  structures  for  instrumentation  and  flight  monitoring 
purposes. 

Because  of  the  importance  and  increasingly  wide  use  of  these  struc¬ 
tures,  major  advances  in  their  analysis  and  design  have  occured  over 
the  past  10  or  15  years.  At  present,  the  Cognac  Platform  off  the 
coast  of  Louisiana  in  the  Gulf  of  Mexico  holds  the  record  as  the 
tallest  offshore  platform.  This  1,265  foot  platform,  which  is  located 
in  1,025  feet  of  water,  is  taller  than  the  Empire  State  Building. 

When  under  full  production  in  1983,  Cognac  is  expected  to  deliver  a 
daily  peak  production  of  50,000  barrels  of  oil  and  150,000,000  cubic 
feet  of  natural  gas. 

Wind,  wave,  and  current  loadings  must  be  considered  in  the  analysis 
and  design  of  fixed  offshore  structures.  Of  these  applied  loads,  wave 
loading  is  by  far  the  most  important  and  least  understood.  The  most 
common  approach  in  determining  wave  loading  is  to  assume  a  regular 
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design  wave  of  height  and  period  representative  of  extreme  storm  con¬ 
ditions.  Calculations  are  then  based  upon  this  wave  for  structure 
analysis  and  design.  A  second  approach  is  to  represent  storm  condi¬ 
tions  using  a  simplified  random  wave  model  and  determine  the  wave  forces 
and  structure  response  in  a  statistical  sense. 

While  both  of  these  methods  are  presently  employed  in  offshore  work, 
the  foundations  on  which  the  analyses  rest  are  still  subject  to  some 
questions.  For  example,  both  methods  employ  a  simplified  wave  force 
formulation  known  as  the  Morison  equation  which  assumes  simply  that 
forces  due  to  water  velocity  and  acceleration  are  additive  at  any 
instant  in  time.  Other  questions  such  as  its  applicability  to  various 
structural  member  orientations  as  well  as  the  proper  choice  of  hydro- 
dynamic  coefficients  indicate  the  uncertainties  still  associated  with 
offshore  structure  design. 

As  a  result,  extensive  theoretical  analysis  of  offshore  structures 
has  occured  over  the  past  few  years,  yet  the  associated  empirical  data 
has  been  limited.  Work  has  been  done  on  wave  forces  on  individual  mem¬ 
bers  as  well  as  limited  work  on  small  scale  structures  in  the  5  to  6 
foot  height  range,  but  facilities  are  generally  unavailable  for  subject¬ 
ing  larger  models  to  both  regular  and  irregular  waves  of  the  proper 
magnitude  and  period  representative  of  storm  conditions.  Limited  results 
have  also  been  reported  on  field  measurements  of  actual  offshore  plat¬ 
forms.  Although  these  tests  are  realistic  for  the  given  sea  state, 
they  lack  the  control  of  waves  offered  by  wave  generation  in  the  labora¬ 
tory. 


The  present  study  was  undertaken  to  investigate  current  techniques 
for  determination  of  wave  force  and  structure  response  under  controlled 
laboratory  conditions.  The  study  is  divided  into  3  areas: 

-  wave  forces  using  regular  waves; 

-  wave  forces  using  irregular  waves; 

-  structural  analysis  using  both  deck  deflections  and 
frequency  response. 


Experimental  Objective 

Experimentally,  the  objective  of  this  study  was  to  measure  the  base 
shear,  lift  force,  overturning  moment,  and  deck  deflection  on  an  off¬ 
shore  platform  test  structure  using  both  regular  and  irregular  waves. 
These  values  are  representative  of  those  required  for  actual  structure 
design  and  analysis  and  were  used  for  comparison  with  associated  theo¬ 
retical  results. 

Testing  was  accomplished  using  the  U.S.  Naval  Academy  Towing  Tank, 
which  measures  380  feet  in  length,  26  feet  in  width,  and  has  a  nominal 
water  depth  of  16  feet.  The  wave  generation  system  is  capable  of  gen¬ 
erating  regular  waves  ranging  from  0.1  Hz  to  1.2  Hz  and  wave  heights  of 
up  to  3  feet.  The  system  is  also  capable  of  generating  periodic  random 
wave  spectra  of  user-defined  specifications.  All  measurements  were  made 
on  a  real  time  basis,  converted  to  digital  format,  and  placed  on  mag¬ 
netic  tape. 


*  * 
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Test  Structure 

The  model  was  designed  and  built  using  thin-walled  6061-T6  aluminum 
tubing  and  is  shown  in  Figures  1  and  2.  The  main  support  piles  con¬ 
sisted  of  2.0  inch  diameter  tubing  of  0.049  inch  wall  thickness  and 
all  cross  bracing  was  of  1.0  inch  tubing  of  0.035  inch  wall  thickness. 
The  design  of  the  test  structure  was  governed  by  3  restrictions: 

-  the  design  was  to  contain  representative  member  bracing  of 
actual  offshore  platforms  while  remaining  simple  enough  to  be  analyzed 
within  the  scope  of  this  project; 

-  member  bracing  was  to  leave  the  structure  flexible  enough 
for  desired  frequency  response  and  deck  deflection  measurements  under 
the  given  wave  loadings; 

-  the  fundamental  frequency  of  the  structure  was  to  be  altered 
by  use  of  added  deck  weight  to  bring  the  structure  response  into  the 
range  for  dynamic  amplification. 

The  structure  was  fixed  to  a  base  assembly  which  allowed  360  degrees 
rotation  of  the  model. 

The  instrumentation  utilized  included  12  4-inch  variable-reluctance 
modular  force  gauges  for  three  dimensional  force  determination,  4  high 
frequency  sonic  transducers  for  deck  displacement  measurements,  and  2 
resistance-type  wave  probes  for  wave  measurements.  A  60-second  computer 
data  acquisition  using  a  sampling  rate  of  51.2  Hz  was  done  on  the  Hydro¬ 
mechanics  Laboratory  Digital  PDP/11  computer  system  for  each  test  run. 

On  each  run  17  data  channels  were  acquired,  and  a  representative  10- 
second  sample  for  regular  waves  and  20-second  sample  for  irregular  waves 
was  truncated  from  all  data  files  for  placement  on  magnetic  tape.  Addi¬ 
tional  test  structure  specifications  and  instrumentation  information  is 
given  in  Appendix  A  and  sample  experimental  data  is  given  in  Appendix  B. 


FIGURE  1  Above  Water  View 


SECTION  I :  REGULAR  WAVE  STUDY 


Numerous  regular  wave  theories  have  been  proposed  ranging  from  a 
simple  linear  theory  developed  by  Airy  in  1845  to  more  complex  non¬ 
linear  theories  such  as  Stokes'  5th  Order  or  Dean's  Stream  Function 
which  have  been  developed  over  the  past  20  years.  Due  to  the  dynamic 
characteristics  of  wave  properties  as  a  function  of  water  depth,  each 
of  these  theories  approximates  the  solution  in  series  form  and  hence 
has  a  region  where  it  gives  best  results.  Thus,  it  is  important  to 
determine  the  type  of  waves  expected  to  act  on  an  offshore  structure 
so  the  proper  wave  theory  may  be  chosen. 

Figure  3  shows  a  regular  wave  moving  at  a  celerity  C  in  water  of 
depth  h.  The  wave  height  H  and  length  L  are  as  shown  and  the  ratio  of 


h 


t 

— > 
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FIGURE  3  Definition  Sketch  for  Regular  Waves 


H/L  is  defined  as  the  wave  slope.  The  wave  period  T  is  related  to  the 
wave  length  and  celerity  by  C=L/T.  The  water  particle  velocity  compo¬ 
nents  at  any  instant  are  u  and  w  and  the  water  surface  displacement 
from  the  still-water  level  (SWL)  at  any  point  is  n. 

The  arrows  at  the  wave  crest,  trough,  and  still -water  lines  indi¬ 
cate  the  directions  of  water  particle  motion  at  the  surface.  These 
motions  cause  a  water  particle  to  move  in  a  clockwise  orbit  as  the 
wave  travels  from  left  to  right.  The  water  particle  velocities  and 
orbit  sizes  decrease  with  increasing  depth  as  shown  in  Figure  4c.  The 
depth  at  which  the  velocity  is  negligible  is  L/2.  A  deep  water  wave 
is  defined  when  h/L^l/2,  transitional  wave  when  l/2<h/L<l/20  and  a 
shallow  water  wave  when  0<h/L<l/20.  As  a  deep  water  wave  travels  toward 
shore,  its  properties  are  altered  as  are  its  particle  orbits  as  shown 
in  Figures  4a  and  4b. 


FIGURE  4  Variation  of  water  particle  orbits  with  depth 
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Linear  Theory 

Linear  wave  theory  was  developed  by  G.  B.  Airy  and  is  significant 
because  of  its  ease  of  application  as  well  as  its  large  range  of  appli¬ 
cability.  Linear  theory  assumes  a  sinusoidal  wave  form  whose  ampli¬ 
tude  is  small  in  comparison  with  its  wave  length  and  the  water  depth. 
Referring  to  Figure  3,  the  free  surface  displacement  from  the  still- 
water  level  is  described  by: 

n  =  H/2  cos(kx-wt)  (1) 

and  the  horizontal  and  vertical  velocity  components  (u,w)  of  the  water 
motion  at  position  (x,z)  and  time  t  are  expressible  as: 

«  '  T  52f®i  <2> 

“  *  f  ^ &W1  <3> 

The  associated  acceleration  components  are  given  by  differentiation 
with  respect  to  time  of  equations  (2)  and  (3).  The  wave  number  k  and 
period  T  are  given  in  terms  of  the  wave  angular  frequency  a>,  and  the 
wave  length  L  by: 

(4) 


T  =  —  ,  k  • 
oi  L 


The  wave  length  L  is  expressible  in  terms  of  the  wave  period  and  water 
depth  as: 

-T?  0_L 

(5) 


L  =  Sli  tanh  L 


2nh 


p  L 

In  deep  water,  however,  tanh  approaches  unity  and  the  wave  length 
can  then  be  expressed  as: 

L  ■  p- 


(6) 
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Although  linear  theory  assumes  a  simple  sinusoidal  wave  form  of  small 
amplitude,  the  agreement  between  theory  and  measured  values  in  deep 
water  is  very  good.  Wiegel  (1964)  showed  that  linear  theory  predicts 
the  velocities  within  the  range  of  experimental  error  for  values  of 
h/L  as  low  as  0.10  for  waves  of  appreciable  steepness.  For  a  more 
complete  description  of  linear  theory,  see  McCormick  (1973)  or  Ippen 
(1966). 

Stokes'  2nd  Order  Theory 

As  the  steepness  of  the  wave  increases,  nonlinearities  become  sig¬ 
nificant  in  the  surface  displacement  and  in  the  velocity  and  accelera¬ 
tion  profiles.  To  account  for  these  nonlinearities,  Stokes  (1880) 
developed  a  wave  theory  in  series  form  and  calculated  second  and  third 
order  corrections  to  the  Airy  theory.  In  deep  water  the  limited  cor¬ 
rections  employed  by  Stokes  are  generally  applicable  but  De  (1955) 
showed  that  Stokes'  theory  should  not  be  used  for  depth  to  wave  length 
ratios  less  than  about  0.125.  A  more  complete  discussion  of  Stokes' 

2nd  order  theory  is  presented  by  Wiegel  (1964)  and  a  brief  summary  of 
the  applicable  equations  for  the  2nd  order  theory  is  given  below. 

In  expanding  to  2nd  order,  the  relation  between  wave  length  and 
period  is  the  same  as  for  the  linear  theory,  namely  equation  (5).  The 
expression  for  the  free  surface  is: 

n  =  i  cos  (kx-ot)  ♦  £-  cosh  fc.h(Z  »  cosh  2kh)  cos  2(kx.„t)  (7) 

2  8L  sinh3  kh 

which  indicates  that  the  mean  level  of  the  wave  surface  lies  above  the 
still-water  level  (i.e.,  the  wave  peaks  more  than  it  troughs).  This 
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is  consistent  with  experimental  data  which  indicate  the  crests  are 

steeper  and  the  troughs  are  flatter  than  those  predicted  by  linear 
theory. 

The  components  of  water  particle  velocities  at  any  point  (x,z) 
and  time  t  in  the  water  are  given  by: 

u  =  S  ".S". *.(.!»»).  cos  (kx-»t)  *  i  (*<f)  gshJMiaii  cos  2(k)t-„t)  (8) 

1  sinh  kh  H  1  L  sinh  kh 

w  =  ^  s,„  (kx.„t)  +  |  (*(*  ^(h*z)  sin  2(kx-„t)  (9) 

sinh  kh  sinh  kh 

and  differentiation  with  respect  to  time  yields  the  water  particle 

accelerations  at  any  point  (x,z). 

Experimental  Waves 

The  regular  waves  chosen  for  this  project  were  of  constant  1/20 
slope  and  were  of  frequencies  from  0.3  Hz  to  1.0  Hz  in  0.1  Hz  incre¬ 
ments.  This  corresponds  to  a  range  of  deep  water  wave  lengths  from 
5.1  feet  to  56.9  feet.  Sorensen  (1978)  shows  that  no  appreciable 
error  will  be  incurred  if  deep  water  is  assumed  for  h/L  ratios  as  low 
as  0.3.  Hence,  with  the  lowest  h/L  ratio  being  0.28,  deep  water  was 
assumed  for  all  theoretical  analyses.  The  choice  of  the  wave  slope 
of  1/20  was  made  in  an  attempt  to  model  storm  waves.  Although  storm 
waves  may  reach  slppe  values  as  high  as  1/10  in  the  Gulf  of  Mexico,  the 
practical  limit  of  the  wave  generation  system,  keeping  wave  distortion 
negligible,  was  a  slope  of  1/20. 


As  shown  in  Figure  5,  Stokes'  2nd  order  theory  is  appropriate  for 
constant  1/20  slope  waves  with  periods  ranging  from  1  to  3.33  seconds 


FIGURE  5  Regions  of  Validity  for  Various  Wave  Theories 

(from  Shore  Protection  Manual ,  U.S.  Army  Corps 
of  Engineers! 
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and  was  therefore  chosen  for  this  study.  Linear  theory  was  used  in 
preliminary  work  and  gave  almost  identical  results. 

Wave  Forces  on  Fixed  Bodies 

The  force  that  a  wave  imparts  on  a  fixed  vertical  pile  was  studied 

by  Mori  son,  et  al.  (1950)  who  proposed  a  semi -empirical  formula  for 

estimating  the  force  per  unit  length  f  in  the  direction  of  the  waves. 

This  has  become  known  as  the  Mori  son  equation: 

2 

f  =  l  pCDD|u|u  +  pCm^-|£  (10) 

where: 

p  =  Density  of  the  fluid 
Co  =  Drag  coefficient 
D  =  Diameter  of  the  pile 
u  =  Water  particle  velocity 
=  Inertia  coefficient 

Jin 

s  Water  particle  acceleration 

This  equation  applies  for  wave  forces  on  cylindrical  surface  piercing 
piles  whose  diameters  are  small  in  comparison  with  the  wave  length. 

The  first  term  in  the  Mori  son  equation  is  a  nonlinear  drag  term 
as  it  is  proportional  to  the  square  of  the  water  particle  velocity. 

The  second  term,  being  proportional  to  the  water  particle  acceleration, 
is  an  inertia  term  which  is  90  degrees  out  of  phase  with  the  drag  term. 
The  dominance  of  either  of  these  terms  is  determined  by  the  wave  charac¬ 
teristics.  If  the  wave  height  H  is  equal  to  or  less  than  the  pile  dia¬ 
meter,  then  the  drag  force  is  negligible.  As  the  wave  height  is 
increased,  however,  the  drag  force  becomes  increasingly  important  to  a 


point  where  the  inertia  term  becomes  negligible.  For  the  waves  used 
in  this  project,  the  drag  force  ranged  from  approximately  30  percent 
to  150  percent  of  the  inertia  force.  This  is  within  the  range  experi¬ 
enced  by  full-scale  offshore  structures  during  storm  conditions  and 
thus  provides  a  good  experimental  basis  for  the  study  of  wave  forces. 

The  values  for  and  CM  to  be  used  in  the  Morison  equation  have 
been  studied  by  Keulegan  and  Carpenter  (1958),  Sarpkaya  (1975), 
Chakrabarti  (1975)  and  others.  Because  waves  are  of  an  oscillatory 
nature,  the  steady  flow  values  of  the  hydrodynamic  coefficients  cannot 
be  used.  The  American  Petroleum  Institute  (1977),  which  summarizes 
current  engineering  practices  for  the  offshore  industry  in  the  United 
States,  recommends  using  a  drag  coefficient  for  circular  cylinders  in 
the  range  from  0.6  to  1.0  and  an  inertia  coefficient  in  the  range  from 
1.5  to  2.0. 

The  values  for  velocity  and  acceleration  for  the  Morison  equation 
are  chosen  from  an  appropriate  wave  theory  as  discussed  in  the  previous 
section.  Because  these  values  decay  with  depth,  the  wave  force  distri¬ 
bution  also  decays  with  depth  as  shown  in  Figure  6,  and  the  total  wave 
force  at  any  instant  can  be  expressed  by  the  integral  equation: 

F  =  JZ 2  f(z)  dz  (11) 

Z1 

The  total  moment  M  about  the  seafloor  (z=-h)  exerted  on  that  part 
of  the  pile  may  similarly  be  given  by: 

M  =  J  2  (z+h)f(z)  dz 


(12) 
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where  z-j  and  z ^  are  the  depths  over  which  the  integration  is  to  be 
performed. 

Wade  and  Dwyer  (1975)  reviewed  4  techniques  presently  used  in  the 
offshore  industry  for  using  Morison's  equation  to  determine  the  wave 
force  on  an  arbitrarily  oriented  cylinder.  It  was  shown  that  these 
techniques  could  lead  to  as  much  as  a  22  percent  difference  in  com¬ 
puted  base  shear,  lift  force,  and  overturning  moment.  The  least 
conservative  of  these  methods  was  presented  by  Chakrabarti,  et  al., 
(1975)  of  Chicago  Bridge  and  Iron  Company.  This  form  was  adopted  for 
use  in  the  present  work. 

The  above  method  involves  resolving  the  water  velocity  and  accelera¬ 
tion  into  components  normal  and  tangential  to  the  cylinder  axis.  The 
normal  component  of  the  wave  force  is  then  calculated  using  the  Mori  son 
equation  which  can,  in  turn,  be  resolved  into  a  three  dimensional  vector 


■>  C 


FIGURE  6  Wave  Force  Decay  With  Depth 


for  engineering  application.  The  tangential  force  components  are 
ignored  for  engineering  calculations  because  the  coefficient  of  drag 
due  to  skin  friction  was  shown  by  Hoerner  (1965)  to  be  approximately 
30  to  120  times  smaller  than  the  drag  coefficient  for  transverse  flow. 

Because  of  the  importance  of  these  equations  to  the  present  study, 
they  are  summarized  below.  Referring  to  Figure  7  and  assuming  a  wave 
propagating  in  the  +x  direction,  the  resulting  water  motion  will  have 
horizontal  and  vertical  velocities  u  and  w  and  corresponding  accelera¬ 
tions  of  ax  and  az-  With  polar  coordinates  6  and  <j>  defining  the 
orientation  of  the  cylinder  axis,  the  magnitude  of  the  water  velocity 
normal  to  the  cylinder  axis  is  given  by: 


FIGURE  7  Member  Axis  Orientation 


and  its  components  in  the  x,  y,  and  z  directions  are  given,  respectively, 
by: 


where: 


un  =  u  -  cx(cxu  *  czw> 


v„  =  -cy  (txu  +  c2w) 


w  =  w  -  c(c  u  +  c,w) 
n  z  x  z 


(14) 


cx  =  si n<j>  cose 

cy  =  sin<)>  sine  (15) 

Cz  =  COS(f> 

The  components  of  the  normal  water  acceleration  in  the  x,  y,  and 
z  directions  are  denoted  by  a  ,  a  ,  and  a  and  similarly  given  by 

1 1 A  1 1  Jr  II*. 

differentiations  of  the  components  of  equation  (14). 

With  these  relations,  the  force  per  unit  of  cylinder  length  acting 
in  the  x,  y,  and  z  directions  is  given,  respectively,  by  the  generalized 
Mori  son  equations: 

fx  =  ?p  CDDvUn  +  pCMITianx 

fy  =  jpCDDwn  +  PCM^ny  (16) 

fz  =  CDDvwn  +  pCM^-anz 

For  a  small  member  such  as  a  horizontal  cross  bracing  on  a  struc¬ 
ture,  where  the  water  motion  does  not  vary  appreciably  over  the  member. 
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average  values  of  u,  w,  a  ,  and  a  may  be  used  in  the  above  equations 

X  z 

and  the  total  force  on  the  member  described  simply  by: 


F  =  f  l 
x  x 


F  =  f  l 

y  y 


fz  =  y 


(17) 


where  t  denotes  the  length  of  the  member.  In  the  more  general  case 
where  the  water  velocities  and  accelerations  do  vary  appreciably  over 
the  length  of  the  member,  the  total  force  must  be  calculated  from  the 
relations: 


F  =  ff  ds  F  =  f  f  ds  F  =  /  f  ds  (18) 

X  x  y  's  y  z  JS  Z 

where  s  denotes  the  distance  along  the  member  axis,  and  the  limits  on 
the  integrals  are  chosen  so  as  to  include  all  of  the  pile  on  which  the 
wave  force  acts.  The  total  moments  about  the  seafloor  may  also  be 
calculated  using  Chakrabarti ' s  formulation  similarly  adapted  to 
equation  (12). 

Computerized  Wave  Force  Deternr nation 

For  deep  water  waves,  such  as  those  used  in  this  study,  where  the 
velocity  profile  is  changing  rapidly  with  depth,  it  is  important  that 
the  above  equations  be  evaluated  using  a  numerical  integration  technique. 
For  the  case  where  an  offshore  platform  is  in  sufficiently  shallow  water 
so  that  the  design  wave  may  be  considered  a  shallow  water  wave,  the 
analysis  may  be  simplified  to  the  more  approximate  approach  of  equation 
(17)  without  incurring  large  errors.  Using  this  approach  on  the  test 
structure  of  the  present  study,  however,  yielded  as  much  as  100  percent 
difference  from  the  more  exact  numerical  integration  technique  ulti¬ 
mately  employed. 


I 
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A  computer  based  three  dimensional  mathematical  model  was  developed 
using  Chakrabarti ' s  formulation  for  arbitrarily  oriented  members  and 
Stokes'  2nd  order  wave  theory.  Computer  tests  were  done  to  deter¬ 
mine  the  maximum  step  size  which  could  be  used  for  a  given  wave  to 
insure  the  numerical  integration  scheme  remained  within  2.0  percent 
of  the  exact  integrated  value.  Input  parameters  included: 

f  =  Wave  frequency  (Hz) 

H  =  Wave  height  (ft) 
a  =  Tower  orientation  (degrees) 

CM  =  Drag  coefficient 
CD  =  Inertia  coefficient 
h  =  Still  water  level  (ft) 

The  program  considered  phasing  distances  of  incremental  member 
length  as  a  function  of  tower  orientation.  Similar  consideration  was 
given  to  moment-arm  distances  as  well  as  incremental  member  orienta¬ 
tion  and  depth.  Extensive  'debugging'  was  done  to  validate  this  model 
which  is  the  basis  for  all  regular  wave  force  data  presented  in  this 
report.  Step  sizes  ranged  from  0.50  inches  to  3.0  inches  and  computer 
run  time  required  between  80  and  160  seconds  per  run.  Each  run  calcu¬ 
lated  the  wave  force  over  1  wave  period. 

The  input  parameters  were  based  on  the  individual  experimental  runs. 
Wave  heights  were  within  +5.0  percent  of  those  required  for  waves  of 
1/20  slope.  One  exception  to  this  was  the  wave  height  for  the  0.3  Hz 
waves  which  were  generally  15  percent  below  the  wave  height  required 
for  1/20  slope.  The  still-water  level  which  changed  3.60  inches  (1.9 
percent)  between  the  2  phases  of  testing  was  determined  to  theoretically 
cause  as  much  as  a  6  percent  error  in  computed  base  shear,  lift  force, 
and  overturning  moment.  In  the  borderline  case  of  deep  water,  this  is 


not  surprising  due  to  the  sensitivity  of  the  wave  equations.  Tank  tem¬ 
perature  variation  was  negligible,  therefore  the  density  of  the  water 
was  assumed  constant. 

Individual  experimental  data  points  in  this  report  represent  the 
average  of  the  peak  values  for  each  run.  An  experimental  error  analysis 
is  given  in  Appendix  C. 

Results  and  Discussion 

The  results  of  the  regular  wave  study  are  presented  in  Figures  8 
through  19  and  are  from  experimental  runs  using  the  previously  mentioned 
waves  for  structure  orientations  of  0  degrees,  30  degrees,  and  45  degrees. 
The  results  of  a  rotational  study  are  also  shown  using  0.5  Hz  waves  of 
0.99  feet  average  height. 

A  range  of  theoretical  results  was  calculated  using  the  above  men¬ 
tioned  computer  program  based  on  a  drag  coefficient  range  of  0.6  to  1.0 
and  inertia  coefficient  range  of  1.5  to  2.0.  This  theoretical  band 
defines  the  range  of  the  American  Petroleum  Institute  (API)  recommended 
standards  and  is  indicated  by  the  shaded  region  of  each  plot. 

The  results  of  the  0  degrees  structure  orientation  study  as  shown 
in  Figures  8,  9,  and  10  indicate  that  experimental  values  of  base  shear 
force,  lift  force,  and  overturning  moment  fall  within  the  theoretical 
limits  defined  by  API  recommended  hydrodynamic  coefficients  for  fre¬ 
quencies  ranging  between  0.3  Hz  and  0.7  Hz.  For  the  higher  frequencies 
ranging  from  0.7  Hz  to  1.0  Hz,  however,  the  experimental  values  consis¬ 
tently  lie  above  the  theoretical  band.  This,  may  be  due  to  member 


phasing  relationships  as  the  wave  lengths  become  increasingly  shorter  or 
to  the  extreme  variability  of  the  drag  and  inertia  coefficients  at  low 
values  of  water  particle  velocity  as  reported  by  Sarpkaya  (1976). 
Reviewing  the  data  points  along  an  entire  curve,  it  can  be  seen  that 
the  position  within  the  theoretical  limits  is  not  constant.  When 
reviewed  on  a  run-by-run  basis  for  the  3  figures,  however,  there  exists 
a  consistency  in  the  position  of  this  point  in  relation  to  the  theore¬ 
tical  limits.  This  may  be  considered  the  result  of  uniform  experimental 
error  for  the  3  measurements  or  more  likely,  the  result  of  variation  of 
actual  drag  and  inertia  coefficients. 

Essentially  identical  results  are  indicated  by  Figures  11  through  16 
which  present  the  results  of  regular  wave  studies  for  structure  orienta¬ 
tions  of  30  and  45  degrees.  Not  only  do  the  experimental  values  follow 
the  same  general  curve,  but  these  points  for  a  given  frequency  can  be 
found  to  lie  in  very  similar  positions  within  the  theoretical  limits. 

This  is  shown  more  clearly  by  plotting  the  results  for  the  rotation 
study  as  shown  in  Figures  17,  18,  and  19.  The  theoretical  values  were 
calculated  over  the  rotation  range  using  a  0.5  Hz  wave  of  0.99  feet  in 
height.  Actual  corresponding  wave  heights  ranged  from  0.97  feet  to 
1.00  feet,  and  the  experimental  data  points  may  be  considered  constant 
within  experimental  error.  The  slight  variation  of  theoretical  values 
from  0  to  45  degrees  (1.8  percent  for  base  shear  force,  10  percent  for 
lift  force,  and  1.8  percent  for  overturning  moment  may  be  due  to  the 
result  of  a  slightly  non-uniform  velocity  profile  acting  across  the 
entire  structure  for  the  given  wave.  If  this  is  true,  replacing  the 
0.5  Hz  wave,  which  has  a  corresponding  wave  length  of  20.5  feet  with  a 
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0.4  Hz  wave,  of  wave  length  32.0  feet,  would  be  expected  to  yield  more 
consistent  theoretical  values  over  the  rotation  range  of  the  structure. 
Using  an  identical  wave  of  0.4  Hz  and  1.64  feet  in  height  and  assuming 
a  constant  water  depth,  the  variation  in  computed  values  over  the 
rotation  range  was  found  to  be  0.05  percent  for  base  shear  force,  0.07 
percent  for  lift  force,  and  0.05  percent  for  overturning  moment,  thus 
substantiating  this  explanation. 

The  consistency  of  the  theoretical  results  for  the  rotation  study 
is  strictly  a  consequence  of  Chakrabarti ' s  formulation  for  computing 
wave  forces.  This  would  not  necessarily  be  true  for  other  methods 
such  as  those  based  on  the  presented  area  and  volume  of  a  structural 
member.  Hence,  the  invariance  of  the  experimental  results  with  rota¬ 
tion  tend  to  substantiate  Chakrabarti ' s  approach  to  wave  force  deter¬ 
mination  on  arbitrarily  oriented  structural  members. 

Although  Figures  17,  18,  and  19  indicate  that  proper  hydrodynamic 
coefficients  should  be  chosen  towards  the  lower  limits  of  the  recom¬ 
mended  API  values,  similar  plots  at  0.3  Hz  would  indicate  near  mean 
values  to  be  appropriate.  Other  plots  at  0.8  Hz  or  0.9  Hz  would  indi¬ 
cate  that  values  at  the  upper  limits  or  even  above  the  recommended  API 
values  would  be  appropriate.  Hence,  although  the  API  values  for  the 
hydrodynamic  coefficients  are  generally  adequate  below  0.8  Hz  for  the 
test  structure,  actual  drag  and  inertia  coefficient  values  appear  to 
have  strict  dependence  upon  some  value  associated  with  frequency  of  the 


water  waves. 
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SECTION  II:  IRREGULAR  WAVE  STUDY 


Advancing  the  wave  force  analysis  a  step  further,  irregular  waves 
may  be  used  to  represent  the  random  nature  of  the  oceans.  A  typical 
wave  trace  of  an  ocean  wave  is  shown  in  Figure  20,  and  because  of 
their  random  nature,  the  characteristics  of  these  waves  are  most 
easily  handled  by  statistics. 

The  simplest  and  most  widely-used  description  of  irregular  waves 
neglects  finite-amplitude  effects  and  uses  linear  theory.  It  assumes 
that  the  surface  deflection  n  from  the  still-water  level  at  any  point 
can  be  represented  by  the  random  sum  of  waves  of  variable  amplitude 
and  random  phase  such  that: 


A„cos(k„x  ~  <*)„t  +  ej 
n  n  n  n 


where: 

U 

A  =  Wave  amplitude  equal  j 
k  =  Wave  number 
u  =  Angular  frequency 
e  =  Phase  angle 

with  the  subscript  indicating  the  value  of  the  nth  wave. 


(19) 


FIGURE  20  Irregular  Wave  Trace 
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The  phase  angle  en  in  equation  (19)  is  assumed  to  be  distributed 
uniformly  over  the  interval  0  to  2*  so  that  any  one  value  is  as  likely 
to  occur  as  any  other.  The  amplitude  An  associated  with  wave  frequen¬ 
cies  between  o>  and  <d  +Aw  is  further  assumed  expressible  as: 
n  n  r 

l A2  =  S  Aw  (20) 

c  r. 

where  is  the  energy  or  amplitude  spectrum  and  represents  the  wave 
energy  distribution  over  the  frequency  interval. 

For  each  wave  component  of  equation  (19),  the  associated  horizontal 
and  vertical  velocities  of  the  water  may  be  determined  from  linear 
theory.  Assuming  the  waves  are  all  propagating  in  the  same  direction, 
the  horizontal  water  particle  velocity  u  and  the  vertical  water  particle 
velocity  w  are  given  at  any  position  x,z  and  time  t  by: 


cosh  kn(h+z) 


u  =  >  U)  A 

- —  n  n 

M 

sinh  kRh 

COS(k  x  -  art  +  e) 
n  n  n 

(21) 

w  =  /_  A 

^ —  n  n 

sinh  kn(h+z) 

sin(knx  -  »nt  +  en) 

(22) 

sinh  k  h 

where  x  is  the  phasing  distance,  z  is  measured  positive  upward  from  the 

still-water  level  and  k  is  related  to  for  deep  water  by: 

n  n 

“n  '  Skn  <23> 

Numerous  theoretical  approaches  have  been  proposed  to  describe  the 
energy  spectrum  of  random  sea  states  based  on  a  combination  of  theore¬ 
tical  reasoning  and  empirical  data.  The  spectrum  used  in  this  report 
was  proposed  by  Bretschneider  in  1959.  This  is  a  2  parameter  spectrum 


based  on  H  and  f  which  may  be  written  as: 
s  m  J 


i  A4 

A  Vr  / 


S  H2e4 

n  4  f5  s 


(24) 


where  Hg  is  the  significant  wave  height,  defined  as  the  average  of  the 
highest  1/3  waves,  and  f  is  the  frequency  of  peak  energy. 


Experimental  Wave  Generation 

Individual  wave  height  components  H(f)  must  b«r  defined  at  descrete 
intervals  of  frequency  so  the  spectrum  may  be  developed  by  computer 
wave  generation.  The  values  of  H(f)  are  a  function  of  the  repeat  period 
which  is  simply  defined  as  the  time  interval  before  the  spectrum  begins 
to  repeat  itself.  The  relation  is: 

1 

H(f)  =  [SS^Af]2  (25) 

where  Af  is  the  reciprocal  of  the  repeat  period.  The  Bretschneider 
spectrum  was  chosen  for  this  project  because  both  the  significant  wave 
height  and  frequency  of  peak  energy  may  be  specified. 

Five  Bretschneider  spectra  were  defined  on  the  computer  with  a  peak 
energy  frequency  of  0.4  Hz  and  a  significant  wave  height  range  of  0.5 
feet  to  1.50  feet  in  0.25  feet  increments.  Individual  wave  components 
were  defined  between  0.3  Hz  and  1.0  Hz  in  0.05  Hz  increments.  Criteria 
for  choosing  these  values  included  keeping  the  peak  energy  frequency 
high  enough  to  assume  deep  water  while  low  enough  to  accomodate  signi¬ 
ficant  wave  heights  up  to  the  1.5  foot  range.  A  plot  of  the  computer 
defined  spectra  is  shown  in  Figure  21. 


Cl 


is  the 


in  which  the  term  | u|  u  is  replaced  by  J 8/tt'  u _ u,  where  u _ 

v  rms  rms 

root-mean  square  of  the  horizontal  water  particle  velocity.  Borgman 
showed  that  this  linearization  represents  the  best  linear  estimate 
in  a  least-square  sense. 

The  restrictions  inherent  in  this  approach  include: 

-  only  vertical  piles  may  be  used  in  modeling  the  structure; 

-  forces  are  only  considered  below  the  still-water  level; 

-  no  vertical  velocity  components  may  appear  in  the  force 


equation. 

This  last  restriction  precludes  the  use  of  a  modified  Chakrabarti  formu¬ 
lation  for  irregular  wave  force  determination  with  Borgman' s  approach. 

Using  this  approximation,  the  Morison  equation  becomes  expressible 
for  irregular  waves  as: 


f  _  v  2  r  ,  y  r  *0  3u 

f  =  —  pCr-Du^u  +  pCu-s-  — r 

n  ir  D  rms  n  M  4  at 


(26) 


where  u  is  dependent  upon  depth  and  is  determined  from  integration 


of  S  by  the  relation: 


U2  ,  £  coshMhWljZ  s  d 

rms  Jo  L  sinh  kh  J  n 


(27) 


Further  using  the  wave  energy  spectrum,  the  force  spectrum  of  the 
structure  may  be  expressed  as: 


sF  =  t;  S 
F  F  n 


(28) 


where  Tf  is  the  transfer  function  relating  Sp  and  and  is  equal  to 
the  maximum  force  experienced  by  the  structure  for  unit  amplitude  waves 


at  specified  frequencies  of  S  . 
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Similarly  from  equation  (28),  the  total  moment  spectrum  SM  of  the 
structure  may  be  expressed  as: 


SM  =  T2  S 
M  M  * 


(29) 


where  is  the  transfer  function  relating  S  with  S 


M' 


Theoretical  Calculations 

The  restriction  of  the  Borgman  method  for  treating  random  wave 
forces  required  the  test  structure  to  be  modeled  using  vertical  cir¬ 
cular  cylinders.  Representative  areas  and  volumes  of  the  horizontal 
and  diagonal  members  were  lumped  at  10  points  along  each  pile.  An 
integration  scheme  for  determining  forces  and  moments  based  on  pre¬ 
sented  area  and  volume  was  developed  which  numerically  integrated 
along  the  length  of  the  piles  and  accounted  for  phasing  relationships. 

As  a  check  on  the  accuracy  of  this  idealized  description,  the  model 
was  used  in  a  regular  wave  computer  study  and  the  maximum  computed  wave 
force  compared  with  that  found  using  the  more  exact  model  employed 
earlier  in  the  regular  wave  study.  The  results  varied  from  -5  percent 
to  -95  percent  for  unit  amplitude  waves  and  frequencies  ranging  from 
0.3  Hz  to  0.9  Hz.  It  should  be  recalled  that  the  total  force  of  the 
structure  using  a  statistical  analysis  depends  upon  u  ,  which  is 
dependent  upon  the  entire  bandwidth. 

To  improve  upon  the  idealized  structural  description,  a  correction 
procedure  was  developed  which  made  the  error  of  the  maximum  structure 
force  for  unit  amplitude  waves  uniformly  small  over  the  bandwidth  while 
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remaining  within  restrictions  imposed  by  the  Borgman  analysis.  The 
computer  program  SMEARTEK  determined,  by  iteration,  the  necessary 
correction  values  of  lumped  areas  and  volumes  based  on  velocity  pro¬ 
file  and  phasing  considerations.  These  values  were  then  implemented 
into  the  previously  mentioned  regular  wave  study  and  the  results 
indicated  the  simplified  model  now  differed  only  by  a  maximum  -18 
percent  for  1  point  with  all  other  values  within  ±5  percent.  The 
results  are  shown  in  Figure  22.  This  'smeared'  model,  when  used  with 
sufficiently  large  irregular  wave  height  components  to  induce  drag 
force  dependency,  provided  a  good  basis  for  experimentally  determining 
the  applicability  of  Borgman 's  analysis. 

Another  application  of  the  program  SMEARTEK  was  to  determine 
correcting  values  for  regular  waves  of  1/20  slope  over  the  frequency 
bandwidth.  Using  these  values,  the  total  base  shear  force  and  over¬ 
turning  moment  were  determined  within  5  percent  of  those  determined 
by  Chakrabarti ' s  formulation.  For  this  application,  the  restriction 
to  linear  wave  theory  and  force  calculations  only  below  the  still- 
water  level  was  no  longer  necessary,  and  values  were  obtained  using 
only  5  percent  to  10  percent  of  the  total  computer  run  time  required 
for  the  more  exact  analysis.  The  results  of  the  SMEARTEK  correction 
for  1/20  slope  waves  are  shown  in  Figure  23. 

An  additional  computer  program  was  developed  to  determine  the 
force  and  overturning  moment  spectral  densities  at  specified  frequen¬ 
cies  using  the  experimentally  determined  wave  spectral  density.  This 
program  was  run  for  both  corrected  and  uncorrected  wave  force  models 
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using  the  Borgman  analysis.  Determination  of  at  each  incremental 
value  of  depth  was  determined  by  integration  of  the  experimentally 
determined  wave  spectral  density  of  Sn  as  given  by  equation  (27)  using  the 
trapazoidal  rule  with  a  Af  of  0.05  Hz.  The  transfer-function  calcula¬ 
tion  considered  regular  waves  of  unit  amplitude  and  horizontal  velo¬ 
cities  u  and  u  consistent  with  the  Borgman  approximation.  The 
actual  value  of  the  transfer  function  was  simply  the  maximum  value  of 
force  or  moment  on  the  idealized  structure  for  individual  frequency 
components. 

Results  and  Discussion 

A  computer-defined  Bretschneider  spectrum  with  a  significant  wave 
height  of  1.00  feet  was  developed  in  the  Naval  Academy's  380  foot 
Towing  Tank.  The  experimentally  determined  wave  spectrum  had  a 
significant  wave  height  of  1.01  feet  and  is  shown  in  Figure  2d.  The 
experimental  force  spectrum  Sp  is  shown  in  Figure  25  along  with  the 
theoretically  determined  force  spectrum  range  using  API  recommended 
drag  and  inertia  coefficients.  The  corresponding  overturning  moment 
spectrum  SM  is  shown  in  Figure  26.  These  figures  indicate  theoretical 
values  using  the  uncorrected  lumped  area  and  volume  technique  for 
cross  bracing.  The  same  experimental  values  are  shown  in  Figures  27 
and  28  with  the  theoretical  range  indicating  values  determined  using 
SMEARTEK  correction  values. 

The  results  using  the  uncorrected  idealized  model  indicate  the 


experimental  values  lie  toward  the  upper  bounds  of  the  theoretical 
range.  The  corrected  structural  model  as  indicated  by  Figures  27  and 
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28  shifted  this  range  to  higher  values.  This  was  expected  due  to  the 

f 

results  of  the  theoretical  regular  wave  tests  run  on  this  model.  The 
degree  to  which  the  range  was  shifted  varied  from  21  percent  at  the 
lower*  frequencies  up  to  59  percent  at  0.9  Hz.  Even  though  the  theo¬ 
retical  values  at  the  upper  frequency  range  increased  more  in  a  per¬ 
centage  sense  as  compared  to  the  lower  frequencies,  the  force  and 
moment  calculations  above  0.6  Hz  were  not  sensitive  enough  to  repro¬ 
duce  the  secondary  'hump'  of  the  wave  spectral  density  and  experimental 
force  spectral  density  which  appeared  between  0.6  Hz  and  0.8  Hz.  This 
corrected  model  appears  to  provide  a  better  theoretical  range  for  the 
experimentally  determined  values. 

Because  irregular  wave  analysis  employs  a  statistical  approach, 
the  analysis  is  subject  to  errors  arising  from  the  following  3  sources: 

-  the  mathematical  approach  suggested  by  Borgman; 

-  the  spectral  determination  method; 

-  the  idealized  model  development  technique. 

The  overall  agreement  between  theoretical  and  experimental  spectral 
values  can  be  considered  within  the  above  stated  restrictions  of  the 
statistical  analysis  as  well  as  any  errors  in  experimental  measurements. 
This  indicates  the  Borgman  method  for  approximating  force  and  moment 
spectral  densities  by  linearizing  the  drag  term  may  be  used  with  a 
degree  of  confidence  in  spectral  density  calculations  of  offshore  plat¬ 
forms. 


FIGURE  22  RESULTS  OF  SMEARTEK  CORRECTION  FOR 
IRREGULAR  WAVE  ANALYSIS 


FREQUENCY  (Hz) 


FIGURE  26  TOTAL  OVERTURNING  MOMENT  SPECTRAL  DENSITY 


FREQUENCY  (Hz) 


FIGURE  28  TOTAL  OVERTURNING  MOMENT  SPECTRAL  DENSITY 
USING  CORRECTED  STRUCTURAL  MODEL 
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SECTION  III:  STRUCTURAL  RESPONSE  STUDY 


Offshore  platforms  consist  of  a  complex  network  of  individual  mem¬ 
bers  designed  to  resist  the  applied  wind,  wave,  and  current  loadings. 
The  response  of  the  entire  structure  as  well  as  individual  member 
response  is  of  primary  importance.  The  analysis  is  carried  out  using 
structural  matrix  methods  which  have  become  practical  for  complex 
structure  analysis  with  the  introduction  of  the  digital  computer. 

From  strength  of  materials,  the  forces  on  a  member  are  proportional 
to  the  member  displacements  and  rotations  in  a  frame  analysis  through 
proportionality  constants  which  express  the  member's  stiffness.  If 
the  member  is  considered  in  two  dimensions,  the  forces  and  moments  at 
each  end  may  be  related  to  the  displacements  and  rotation  of  these 
end  points  by: 


where: 


E  =  Modulus  of  elasticity 
I  =  Moment  of  inertia 
A  =  Member  cross  sectional  area 
l  =  Member  length 

If  the  same  member  is  considered  in  three  dimensions,  a  similar 
relation  to  equation  (30)  may  be  developed  which  includes  axial  and 
transverse  loadings  of  the  member  and  twisting  and  bending  about  2 
principle  axes  of  the  cross  section.  This  necessitates  that  equation 
(30)  be  expanded  to  a  twelfth  order  system  of  equations. 

Considering  the  two  dimensional  case,  equation  (30)  may  be  rewritten 
as : 

{F}  =  [K]  {U}  (31) 

where  the  elements  of  { F } ,  [K],  and  {11}  are  given  by  the  corresponding 
elements  of  equation  (30).  The  matrix  [K]  is  referred  to  as  the  stiff¬ 
ness  matrix  of  the  member. 

To  allow  for  the  specification  of  either  moments  at  the  joints  or 
rotations  of  the  joints,  the  structure  is  usually  modeled  as  a  rigid 
joint  frame.  This  requires  the  joints  to  preserve  the  relative  angular 
spacing  of  the  attached  members  at  the  joints  while  undergoing  loading. 
With  rigid  joints,  all  members  connected  at  the  joint  will  therefore 
suffer  the  same  rotation  at  the  joint  as  the  joint  itself. 

Multi-member  structures  such  as  offshore  platforms  are  analyzed  in 


this  manner  using  properly  constructed  matrices  which  reflect  individual 
member  connectivity,  properties,  and  orientation.  This  method  of 
structural  analysis  has  become  known  as  the  direct  stiffness  method. 


57 


Regardless  of  the  size  of  the  structure,  equation  (31)  always  applies, 
the  only  difference  being  the  size  of  the  matrices.  For  example,  the 
test  structure  shown  in  Figure  29a  in  two  dimensions  requires  a  36th 
order  matrix  to  reflect  2  degrees  of  translation  and  1  degree  of  rota¬ 
tion  at  each  joint.  The  same  structure  in  three  dimensions  as  shown  in 
Figure  29b,  however,  requires  a  144th  order  matrix  to  handle  the 
3  degrees  of  freedom  translation  and  3  degrees  of  freedom  rotation 
wr.ich  is  necessary  to  describe  the  system. 

Applying  the  proper  boundary  conditions  to  insure  a  determinate 
system  and  knowing  the  forces  and  moments  at  each  joint,  equation  (31) 


a.  Two  Dimensions  b.  Three  Dimensions 

FIGURE  29  Computer  Based  Math  Models 
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may  be  solved  for  the  displacements  and  rotations  through  the  relation: 

{U}  =  [K]"1  {F}  (32) 

where  [K]  ^  is  the  inverse  of  the  stiffness  matrix  or  more  commonly 

called  the  flexibility  matrix. 

Determination  of  the  forces  and  moments  on  each  member  is  accom¬ 
plished  through  the  use  of  an  appropriate  wave  and  force  theory  as 
described  in  previous  sections.  As  a  wave  passes  the  structure,  there 
are  2  cases  of  importance  -  the  first  being  when  the  total  force  on 
the  structure  is  at  a  maximum  and  the  second  being  when  the  force  on 
an  individual  member  is  at  a  maximum.  Because  the  overall  structure 
response  as  indicated  by  deck  deflection  is  of  interest  in  this  project, 
only  the  first  case  was  studied. 

When  the  total  structural  force  is  at  a  maximum,  individual  member 
forces  may  be  computed  and  modeled  using  2  techniques.  The  first  tech¬ 
nique  as  shown  in  Figure  30a  calculates  the  force  on  the  member  and 
distributes  this  load  equally  over  the  member.  This  load  may  then  be 
computed  into  equivalent  point  loads  and  moments  as  shown  in  Figure  30b 


a.  distributed  load 
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b.  equivalent  loads 


FIGURE  30  Distributed  Load  Technique 


and  applied  at  the  end  points.  The  second  technique  assumes  a  concen¬ 
trated  load  acting  on  the  member  and  similarly  computes  equivalent  joint 
loads  and  moments  as  shown  in  Figures  31a  and  31b. 

Using  a  digital  computer,  equation  (32)  may  be  solved  for  the  unknown 
displacements  and  rotations  of  each  joint.  These,  in  turn,  may  be  used 
to  determine  the  internal  forces  in  each  member  from  which  the  maximum 
stresses  may  be  found  from  well  established  solid  mechanics  principles. 

For  a  complete  discussion  of  structural  analysis  by  the  direct  stiffness 
method  see  Willems  and  Lucas  (1968). 

Computer  Solution 

The  matrices  described  by  equation  (31)  were  developed  fora  two  dimen¬ 
sional  model  of  the  test  structure  and  placed  on  the  computer.  The  pro¬ 
gram  used  for  the  regular  wave  analysis  was  then  adapted  to  use  the 
concentrated  load  technique  described  above.  The  method  of  distributed 
loads  was  once  again  found  undesirable  due  to  the  deep  water  velocity 
profiles  of  the  experimental  waves.  The  step  size  considerations  remained 
the  same  as  in  the  regular  wave  force  analysis  to  insure  accurate  results. 


P 


a.  concentrated  load  b.  equivalent  loads 


FIGURE  31  Concentrated  Load  Technique 


To  verify  the  concentrated  load  technique  used  in  this  program, 
the  results  were  compared  with  the  regular  wave  force  program.  The 
maximum  percentage  difference  for  all  encountered  waves  was  7.5  per¬ 
cent.  This  error  could  be  attributed  to  the  uncertainties  of  the 
technique  for  diagonal  members  normal  to  the  waves. 

The  equivalent  forces  and  moments  at  each  joint  found  by  this 
method  were  then  used  in  the  previously  constructed  matrices.  The 
matrix  equation  was  solved  on  the  computer  for  individual  joint  dis¬ 
placements  and  rotations. 

Frequency  Analysis 

As  offshore  structures  move  into  deeper  water,  they  tend  to  become 
more  flexible  and  as  a  result  the  dynamic  aspects  of  wave  loadings 
become  more  important.  The  natural  frequencies  of  these  platforms  can 
be  excited  if  the  frequency  of  the  design  wave  or  peak  frequency  of 
the  design  spectrum  comes  within  range  of  the  platform  frequencies. 

The  design  of  the  test  model  required  2  frequency  considerations; 
the  first  was  to  have  the  fundamental  response  high  enough  to  assume 
static  wave  loading  and  the  second  was  to  allow  this  response  to  be 
changed  sufficiently  using  added  deck  weight  to  allow  a  study  of  dynami 
loading  effects. 

Single  Spring/Mass  Model 


Using  the  previously  constructed  matrices  for  the  two  dimensional 
structural  analysis,  a  load  evenly  placed  on  the  deck  of  the  platform 


will  result  in  a  given  deck  displacement  such  that  the  structure's 
effective  spring  constant  K'  may  be  calculated  from  the  relation: 

K'  =  £  (33) 

where  U  is  the  deck  displacement  for  the  load  F.  Assuming  a  1  degree 
of  freedom  spring/mass  system,  the  first-mode  natural  frequency  of  the 
structure  may  be  calculated  from: 


where  M'  can  be  considered  the  effective  mass  of  the  structure  given 
approximately  by: 

M*  =  M  +  ^Mp  +  1mv  (35) 

where: 

M  =  Mass  of  the  deck 

Mp=  Mass  of  the  support  structure 

Mv=  Mass  of  the  displaced  water 

This  method  uses  a  simple  lumped  mass  model  which  approximates  for  the 
mass  of  the  structure  and  displaced  water. 

When  the  above  method  was  applied  to  the  two  dimensional  test  model, 
the  fundamental  frequency  response  was  found  to  be  4.66  Hz  for  a  deck 
weight  of  125  pounds. 


Three  Dimensional  Analysis 


The  complexity  of  a  frequency  analysis  fora  three  dimensional  platform 


has  lent  itself  to  the  development  of  computer  codes  to  handle  user 


defined  platforms.  Modeling  these  structures  as  multi-degrees  of 
freedom  systems  allows  for  a  more  accurate  and  complete  frequency 
analysis. 

The  computer  code  'GIFTS-IV*  (Graphics  oriented  Interactive 
Finite  element  Time  sharing  System)  was  used  to  model  the  test 
structure  as  a  three  dimensional  multi -degrees  of  freedom  system 
for  the  calculation  of  6  modal  frequencies.  'GIFTS'  was  developed 
by  Professor  Hussein  A.  Kamel  and  Michael  W.  McCabe  of  the  Univer¬ 
sity  of  Arizona  as  a  modular  code  for  structural  analysis.  The 
structure  to  be  analyzed  must  first  be  defined  and  may  then  be 
subjected  to  both  static  and  transient  load  analyses  as  well  as  a 
frequency  analysis.  The  test  model  was  developed  in  three  dimen¬ 
sions  as  shown  in  Figure  29b,  and  the  first  mode  frequency  analysis 
as  described  by  equation  (34)  resulted  in  a  response  within  1  per¬ 
cent  of  the  two  dimensional  model.  This  is  not  surprising  due  to 
the  symmetry  of  the  model.  The  front  and  rear  faces  contribute 
very  little  to  the  stiffness  of  the  model  for  a  force  applied,  say, 
to  the  positive  x  axis  for  a  tower  orientation  of  0  degrees. 

For  the  three  dimensional  analysis,  'GIFTS’  lumps  the  masses  of 
the  members  at  each  of  the  24  connecting  nodes,  applies  loads  in  the 
x,  y,  and  z  directions  represented  by  a  power  series,  and  then  uses  a 
subspace  iteration  routine  to  converge  on  the  frequency  modes.  This 
routine  in  three  dimensions  produced  results  which  differed  from  the 
experimental  values  of  the  fundamental  mode  by  150  percent  or  more 


(1.38  Hz  predicted  value  versus  3.61  Hz  experimental  for  125  pounds 
deck  weight).  Using  'GIFTS'  for  the  structure  defined  in  two  dimensions, 
however,  yielded  a  fundamental  response  equal  to  that  obtained  by  an 
equation  (34)  analysis.  The  'GIFTS'  frequency  analysis  in  three  dimen¬ 
sions  indicated,  however,  that  the  first  3  modes  of  the  structure 
would  be  almost  identical.  This  corresponded  to  the  experimentally 
observed  cross-coupling  of  these  modes. 

Dynamic  Response 

Another  structural  consideration  which  is  dependent  upon  both  the 

frequency  response  of  the  structure  as  well  as  the  wave  force  loading 

is  dynamic  response.  Regular  wave  loading  on  an  offshore  structure 

can  be  considered  an  oscillatory  forcing  function  of  some  angular 

frequency  u>.  Extending  the  single  spring/mass  model  of  the  test 

structure  to  include  a  dashpot  to  account  for  damping,  the  amplitude 

of  the  deck  deflection  U  as  a  function  of  —  ,  where  <o  is  the  undamped 

“’n  n 

natural  frequency  of  the  structure,  may  be  expressed  in  nondimensional 
form  as: 

^  =  — - 1 - -  (36) 

Fo  •  /[  1  -  (£)2]2  +  [2C(£)]2 

V  wn  wn 

where 

c  =  Damping  factor 

K  =  Structure  stiffness 

F0  =  Magnitude  of  forcing  function 

I  IK 

This  technique  is  discussed  by  Thompson  (1972)  and  the  value  of  has 

ho 

become  known  as  the  dynamic  amplification  factor. 


64 


The  curve  of  the  dynamic  amplification  factor  for  the  test  structure 

as  a  function  of  —  was  generated  on  the  computer.  Because  c  is  an 
n 

experimentally  determined  parameter,  was  taken  to  be  the  experimental 
value  instead  of  the  theoretical  value  described  by  equation  (34)  for 
the  sake  of  consistency. 

Using  the  concentrated  load  technique  described  for  structural 
analysis  but  lumping  the  forces  at  the  deck  and  seafloor  instead  of 
at  each  joint  will  give  the  necessary  Fq  to  solve  equation  (36)  for 
deck  displacements.  The  value  of  K  in  this  case  is  the  effective 
spring  constant  of  the  structure  determined  from  equation  (33). 

Experimental  Determination  of  the  Fundamental  Response 

The  first  mode  response  was  determined  experimentally  using  a  number 
of  excitation  techniques  to  assure  repeatability.  Accelerometers  were 
mounted  on  the  deck  of  the  structure  as  well  as  the  top  horizontal  cross 
members.  A  'twang'  test  and  impulse  test  were  used  for  first  mode 
response  determination.  In  addition,  the  structure  was  pushed  (on  axis) 
and  then  released.  All  3  methods  gave  identical  results  within  the 
accuracy  of  the  experiment  and  the  average  values,  for  the  structure  in 
water,  along  with  the  three  dimensional  theoretical  response  from  equa¬ 
tion  (34).  are  shown  in  Figure  32. 

In  determining  the  frequency  response  of  the  structure,  the  effect 
of  damping  was  assumed  negligible.  This  may  be  checked  using  the  method 
of  logarithmic  decrement  as  discussed  by  Thompson  (1972).  The  damping 
ratio  c  was  determined  from  an  accelerometer  trace  to  be  0.045,  which 
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indicated  that  the  frequency  damping  of  the  structure  was  In  fact 
negligible  =  0.999^).  The  additional  damping  due  to  the  water 
was  determined  experimentally  to  be  approximately  1.0  percent  which 
did  not  contribute  significantly  to  the  frequency  damping  of  the 
structure. 

Two  cases  of  deck  weight  were  used  during  testing  to  meet  the  2 
criteria  discussed  earlier;  125  pounds  and  1,330  pounds  which  gave 
respective  experimental  first  mode  responses  of  3.61  Hz  and  1.27  Hz. 

Results  and  Discussion 

A  plot  of  the  experimental  frequency  results  (in  water)  as  well 
as  the  results  of  a  three  dimensional  equation  (34)  analysis  were  given 
in  Figure  32.  The  experimental  values  are  shown  to  range  from  16 
percent  to  25  percent  lower  than  the  theoretical  values.  With  the 
damping  negligible,  the  most  probable  explanation  lies  with  the 
method  of  theoretical  analysis.  The  assumption  of  rigid  joints  used 
in  determining  the  effective  spring  constant  may  model  the  structure 
too  'stiffly,'  hence  accounting  for  the  lower  experimental  frequency 
values. 

The  results  of  the  deck  deflection  calculations  and  corresponding 
experimental  values  indicated  a  similar  error  in  the  theoretical 
analysis.  As  shown  in  Figure  33,  the  experimental  deck  deflections 
differed  by  a  factor  of  2.6  up  to  as  much  as  4.0  above  the  theoretical 
values.  The  maximum  possible  dynamic  amplification  for  this  weight 
case  was  4  percent  and  hence  did  not  contribute  significantly  to 
these  factors  of  error.  Therefore,  with  the  previously  established 
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accuracy  of  wave  force  calculations  as  well  as  the  accuracy  of  the  two 
dimensional  versus  three  dimensional  structure  analysis,  the  possibility 
of  error  within  the  assumption  of  rigid  joints  is  further  substan¬ 
tiated. 

The  results  of  the  dynamic  amplification  study  as  shown  in 
Figure  34  indicate  the  experimental  values  based  on  deck  deflection 
are  in  reasonable  agreement  with  the  dynamic  amplification  curve  for 
a  structure  damping  ratio  c  of  0.045.  These  values  were  based  on  the 
ratio  of  the  deck  deflection  for  1/20  slope  waves  with  a  deck  weight 
of  125  pounds  to  those  for  1/20  waves  with  a  deck  weight  of  1,330 
pounds.  The  error  from  the  theoretical  curve  ranged  from  1  percent 
at  0.3  Hz  up  to  12  percent  at  0.7  Hz. 


—  Theoretical 
+  Experimental 


FIGURE  33  DECK  DEFLECTION 


CONCLUSIONS 


Based  on  the  material  presented  in  the  above  three  sections,  the 
following  conclusions  may  be  drawn: 

•  The  regular  wave  study,  which  used  deep  water  constant  slope 
waves  at  various  incident  wave  angles  on  the  structure,  indicated 
that  force  calculations  based  on  consideration  of  the  velocity  and 
acceleration  components  normal  to  the  individual  members  yielded 
uniformly  good  results  when  compared  with  experimental  measurements 
for  various  structure  orientations  relative  to  the  incident  wave. 

•  As  offshore  structures  continue  to  move  into  deeper  water 
where  even  the  longest  design  waves  must  be  considered  deep  water 
waves,  an  approach  such  as  the  one  used  in  this  study  should  be  con¬ 
sidered  which  combines  the  above  formulation  for  wave  force  analysis 
with  an  accurate  numerical  integration  technique. 

•  The  examination  of  the  Borgman  irregular  wave  force  analysis 
indicated  good  agreement  with  experimentally  determined  values.  This 
examination  provided  a  better  check  on  the  linearization  of  the  drag 
term  than  Borgman  himself  provided  due  to  the  significant  wave  height 
used  in  relation  to  the  pile  diameter.  For  a  structure  with  signifi¬ 
cant  cross  bracing,  an  alternative  approach  to  a  simple  lumped  area 
and  volume  approximation  of  non-vertical  members  which  accounts  more 
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accurately  for  their  wave  force  and  moment  contributions,  appears  to 
provide  better  spectral  representation  of  the  total  force  and  moment 
on  the  structure. 

•  The  American  Petroleum  Institute  recommended  range  of  values 
for  the  drag  and  inertia  coefficients,  0.6  to  1.0  and  1.5  to  2.0, 
respectively,  provided  theoretical  values  which  encompassed  experi¬ 
mental  values  for  wave  frequencies  below  0.8  Hz  for  the  test  structure. 
The  regular  wave  study  indicated  a  strict  consistency  in  the  relation 
of  the  experimental  values  to  this  theoretical  range.  This  was  con¬ 
cluded  to  be  the  result  of  the  variation  in  values  of  drag  and  inertia 
coefficients  for  low  values  of  frequency  dependent  parameters  such  as 
those  discussed  by  Sarpkaya  (1976)  and  Keulegan  and  Carpenter  (1958). 

At  the  time  of  this  report,  theoretical  force  and  moment  calculations 
using  values  of  drag  and  inertia  coefficients  based  explicitly  on 
encountered  frequencies  have  not  been  determined. 

•  The  results  of  the  structural  response  study  indicated  that 
the  joints  of  the  test  structure  were  overly  restrained  in  the  computer 
model  which  employed  a  rigid  joint  frame  analysis  for  both  frequency 
and  deck  deflection  calculations. 
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TEST  STRUCTURE  AND  INSTRUMENTATION 


In  addition  to  the  specifications  of  the  test  structure  which  were 
presented  in  the  text  of  the  report,  material  tests  were  conducted  to 
determine  the  best  weld  configuration  and  loss  of  yield  strength  due  to 
the  tungsten  inert  gas  (TIG)  welding  process.  The  results  indicated  a 
butt-butt  filet  weld  provided  the  best  weld  continuity  while  reducing 
the  yield  strength  approximately  50  percent  to  a  value  of  17,000  psi. 
Corrosion  tests  were  also  completed  which  indicated  a  submergence  time 
of  approximately  1  month  would  result  in  negligible  corrosion  of  the 
structure. 

The  instrumentation  package  chosen  included  12  4-inch  Hydronautics 
variable-reluctance  modular  force  gauges  to  measure  forces,  4  Wesmar 
high  frequency  sonic  transducers  to  measure  displacements,  2  MTS 
resistance-type  wave  gauges  to  measure  wave  characteristics,  and  asso¬ 
ciated  equipment. 

Three  force  gauges  were  used  on  each  leg  as  shown  in  Figures  A-l 
and  A-2  to  allow  the  measurement  of  the  force  in  the  x,  y,  and  z  direc¬ 
tions.  Signals  were  processed  through  Hydronautics  Multi-T  signal  condi¬ 
tioning  units  and  filtered  through  Ithaco  filters  at  10  Hz.  A  complex 
design  of  multi-bearing/pivot  connectors  also  shown  in  Figures  A-l  and 
A-2  was  developed  and  built  to  allow  the  attachment  of  the  structure  to 
the  force  gauges.  This  assembly  allowed  negligible  pre-stressing  of  the 


structure/gauge  interface  which  could  result  from  misalignment  of  the 
structure,  base  plate  bending  and/or  thermal  contraction.  This  system 
was  necessary  to  resolve  accurately  the  force  gauge  outputs  into  resul¬ 
tant  forces  and  moments  on  the  structure. 

Deck  deflection  and  rotation  measurements  were  accomplished  using 
high  frequency  sonic  transducers  in  the  configuration  shown  in  Figure 
A-3.  Signals  were  processed  through  Wesmar  signal  conditioning  units 
and  filtered  through  Ithaco  filters  at  10  Hz.  The  deck  deflection  sig¬ 
nals  were  digitally  manipulated  to  achieve  x  and  y  axis  deflection  and 
deck  rotation.  The  high  speed  towing  carriage  was  placed  over  the 
structure  as  shown  in  Figure  1  to  act  as  the  reference  base  for  the 
deck  deflection  measurements  as  well  as  the  relay  point  between  the 
gauge  signals  and  the  computer. 


FIGURE  A-l  Structure  -  Gauge  Interface 


I 


Wave  characteristics  were  recorded  using  a  resistance-type  wave 
probe  aligned  at  the  leading  edge  of  the  structure  (base  level).  Th 
corresponded  to  the  reference  point  used  in  theoretical  calculations 
and  allowed  for  phasing  comparisons. 


FIGURE  A-2  Structure  -  Gauge  Interface 


tank  axis 


top  view 


f IGURE  A-3  Deck  Deflection  Measurement 'Schematic 


APPENDIX  B 

SAMPLE  EXPERIMENTAL  DATA 


81 


SAMPLE  EXPERIMENTAL  DATA 


Seventeen  channels  of  data  were  acquired  by  the  computer  on  each 
test  run.  Of  these,  12  channels  were  force  data,  4  channels  were 
displacement  data,  and  1  channel  was  wave  data.  In  addition  to 
this  primary  data,  up  to  8  channels  were  calculated  by  digital 
manipulation  of  the  primary  data  files.  Values  such  as  deck  rotation 
and  torsional  moment  were  found  either  too  small  to  measure  or  gave 
meaningless  results  and  are  not  presented.  The  computed  values  pre¬ 
sented  in  this  report  are  limited  to  total  base  shear,  lift  force, 
overturning  moment,  and  x-axis  deck  deflection. 

The  sample  data  presented  in  this  appendix  is  from  both  regular 
and  irregular  wave  test  runs. 


Regular  Wave:  Phase  I  Run  21  0.4  Hz  1.56  ft 

Structure  Orientation  -  Zero  Degrees 

Figure  B-l  Wave  0.4  Hz  1.56  ft 

Figure  B-2  X-axis  force  component  on  leg  'A' 

Figure  B-3  Y-axis  force  component  on  leg  'A' 

Figure  B-4  Z-axis  force  component  on  leg  'A1 

Figure  B-5  X-axis  deck  deflection 

Figure  B-6  Total  base  shear  force 

Figure  B-7  Total  lift  force 

Figure  B-8  Total  overturning  moment 

Irregular  Wave  Phase  II  Run  95  fm  =  0.4  Hz  1.01  ft 

Structure  Orientation  -  Zero  Degrees 

Figure  B-9  Wave  fm  =  0.4  Hz  Hs  =  1.01  ft 

Figure  B-10  X-axis  force  component  on  leg  'A' 

Figure  B-l 1  Y-axis  force  component  on  leg  'A' 

Figure  B-l 2  Z-axis  force  component  on  leg  'A* 


Figure  B-13  Total  base  shear  force 

Figure  B-14  Total  overturning  moment 

Figure  B-15  Spectral  density  of  wave 

Figure  B-16  Spectral  density  of  base  shear  force 

Figure  B-17  Spectral  density  of  overturning  moment 


STRUCTURE  OR 


X-AXIS  FORCE  COMPONENT  ON  LEG  'A' 
STRUCTURE  ORIENTATION  -  ZERO  DEGREES 


FIGURE  B-4  Z-AXIS  FORCE  COMPONENT  ON  LEG  'A' 

STRUCTURE  ORIENTATION  -  ZERO  DEGREES 


FIGURE  3-5  X-AXIS  DECK  DEFLECTION 

STRUCTURE  ORIENTATION  -  ZERO  DEGREES 


FIGURE  B-7  TOTAL  LIFT  FORCE 

STRUCTURE  ORIENTATION  -  ZERO  DEGREES 


1.4 


STRUCTURE  ORIENTATION  -  ZERO  DEGREES 
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FIGURE  B-10  X-AXIS  FORCE  COMPONENT  ON  LEG  'A' 

STRUCTURE  ORIENTATION  -  ZERO  DEGREES 
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Y-AXIS  FORCE  COMPONENT  ON  LEG  'A' 
STRUCTURE  ORIENTATION  -  ZERO  DEGREES 


TOTAL  BASE  SHEAR  FORCE 

STRUCTURE  ORIENTATION  -  ZERO  DEGREES 


TOTAL  BASE  SHEAR  SPECTRAL  DENSITY 
STRUCTURE  ORIENTATION  -  ZERO  DEGREES 


TOTAL  OVERTURNING  MOMENT  SPECTRAL  DENSITY 
STRUCTURE  ORIENTATION  -  ZERO  DEGREES 
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EXPERIMENTAL  ERROR  ANALYSIS 


The  purpose  of  this  analysis  is  to  identify  and  attempt  to  quantify 
those  areas  which  effect  the  accuracy  of  the  experimental  results  in 
this  report. 

Four  areas  can  be  identified  which  could  contribute  to  experimental 
error: 

-  Model  construction; 

-  Physical  parameters; 

-  Measuring  equipment/calibration; 

-  Data  repeatability. 

Model  Construction 

The  construction  accuracy  of  the  18  foot  test  structure  was  excel¬ 
lent.  Dimensions  were  specified  in  the  design  within  ±  1/64  inch  with 
the  desired  construction  accuracy  at  ±  1/4  inch.  Measurements  taken 
of  the  test  structure  after  completion  indicated  all  but  3  measurements 
were  within  specified  limits.  The  3  values  outside  this  range  were 
within  +  3/8  inch.  These  values  were  taken  into  consideration  for  the 
computer  math  models. 

Physical  Parameters 

As  discussed  in  the  text  of  the  report,  the  still-water  level  varied 
3.60  Inches.  This  was  determined  to  vary  the  theoretical  calculations 


as  much  as  6  percent  for  a  0.3  Hz  wave  and  as  such  was  included  as  an 
input  parameter  for  all  theoretical  calculations. 


Tank  temperature  varied  7.8  degrees  Fahrenheit,  which  had  a  negli 
gible  effect  on  water  density  and  as  such  was  not  considered  in  theo¬ 
retical  calculations. 

The  effect  of  tank  seiche  on  experimental  accuracy  was  checked  by 
running  a  0.5  Hz  wave  of  0.79  feet  for  2  minutes  in  the  tank.  The 
tank  was  allowed  to  settle  for  15  minutes  and  a  data  acquisition  was 
taken  for  all  channels  which  indicated  a  natural  tank  period  of 
approximately  32  seconds  with  a  residual  H^  of  0.17  inches.  The 
force  blocks  for  all  3  axes  indicated  force  levels  well  within  the 
noise  level  of  the  signals.  Hence,  the  effects  of  tank  seiche  on 
experimental  values  given  in  this  report  were  considered  negligible. 


Measuring  Equipment/Calibration 


The  electronic  equipment  used  in  this  project  was  based  on  a 
voltage  range  of  ±10  volts.  All  measuring  equipment  was  calibrated 
before  testing  and  signals  were  checked  for  mechanical  zero  on  a 
daily  basis.  A  check  on  the  wave  probe  calibration  was  also  done 
on  a  daily  basis. 

Although  there  were  2  types  of  measurements  on  the  structure, 
deck  deflection  and  forces,  electronically  they  were  similar.  The 
following  is  a  brief  description  of  each  piece  of  equipment  In  the 
signal  path.  Specific  values  not  normally  associated  with  the 
general  equipment  are  given  for  the  particular  application  in  this 
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-  Hydronautics  Modular  Force  Gauges 

These  gauges  utilize  variable  reluctance  displacement  trans¬ 
ducers  to  accurately  measure  an  applied  load.  The  gauges  convert 
the  relative  displacement  of  the  parallel  faces  of  the  cube  to  gen¬ 
erate  an  output  voltage  that  is  a  linear  function  of  the  applied 
load.  They  are  essentially  insensitive  to  cross-axis  loading. 

The  gauges  have  a  linearity  of  ±  0.25  percent  of  design  load  (250 
pounds)  and  a  corresponding  hysterisis  of  0.05  percent  of  the  design 
deflection  (.05  inches).  The  actual  resolution  of  the  gauge  was  0.3 
percent  of  the  load. 

-  Wesmar  LM4000  High  Frequency  Sonic  Transducers 

This  system  can  be  used  for  in-air  range  measurements  from  2 
to  7.5  inches  with  a  resolution  better  than  0.5  percent  of  measured 
range.  The  linearity  is  better  than  0.5%  of  the  full  scale  measure¬ 
ment.  The  response  time  for  the  present  study  was  79.5  Hz  with  a 
half  power  beam  width  of  3  degrees. 

-  Ithaco  4120  Series  Dual  Low-pass  Filters 

These  are  precision  active  filters  with  phase  and  amplitude 
tracking  between  filters  with  the  same  setting  within  ±  4  degrees  and 
±  0.2  dB  respectively.  For  this  project,  the  filters  can  be  considered 
to  have  had  infinite  resolution  and  essentially  zero  nonlinearity  in 
the  pass  band. 

-  Hydronautics  Multi-T  Signal  Conditioner  Unit  (SCU) 

These  are  multi -transducer  conditioners  which  incorporate  the 
means  to  balance,  establish  polarity,  excitation,  calibration,  and 
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to  amplify  signals  from  various  transducers. 

When  used  in  the  reluctance  mode,  the  carrier  frequency  is  1000 
Hz  at  5.5  VRMS  with  an  output  linearity  of  ±  0.2  percent  of  full  scale 
voltage. 

-  Preston  32  Channel  A/D  Converter 

For  the  present  study  the  A/D  Converter  had  a  14  bit  magnitude 
resolution  capability  which  translated  to  an  accuracy  of  0.01  percent 
at  ±  (1/2  LSB)  for  10  volts. 

The  overall  result  indicated  that  the  measuring  equipment  could 
resolve  the  measurements  within  the  accuracy  of  the  calibration  devices. 

Force  gauges  were  calibrated  using  precision  weights  accurate  to 
within  +0.5  percent.  One  pound  increments  were  used  up  to  ±  5  pounds 
and  the  calibration  was  extended  to  ±  20  pounds  using  5  pound  weights. 
Maximum  individual  force  values  were  ±  10  pounds  on  x-axis  gauges, 

±  3  pounds  on  y-axis  gauges,  and  ±  30  pounds  on  z-axis  gauges.  The 
sonic  transducers  were  calibrated  using  0.25  inch  aluminum  plates  with 
an  expected  accuracy  of  ±  1.0  percent.  The  sonic  calibrations  were 
checked  on  a  daily  basis.  The  wave  probe  was  also  checked  on  a  daily 
basis  with  typical  accuracy  being  within  1.0  percent  for  a  12  inch  span. 

Data  Repeatability 

The  data  was  checked  for  repeatability  over  the  entire  60  second 
acquisition,  with  a  representative  10  second  sample  being  placed  on 
magnetic  tape.  In  addition,  identical  runs  at  0.3  Hz  with  a  structure 


orientation  of  0  degrees  indicated  a  2.7  percent  difference  in  total 
base  shear  force,  3.1  percent  in  lift  force,  and  2.6  percent  in  total 
overturning  moment. 
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APPENDIX  D 
COMPUTER  PROGRAMS 


COMPUTER  PROGRAMS 


This  project  resulted  in  the  development  of  35  computer  programs, 
31  by  the  author  to  handle  all  engineering  computations  and  4  by 
Mr.  Gene  Morgan  and  Mr.  Lon  Ward  of  the  Applications  Group,  Academic 
Computing  Center,  to  handle  the  transfer  and  digital  manipulation  of 
experimental  data.  In  addition,  7  programs  were  adapted  by  Dr.  John  S. 
Kalme,  Mathematics  Department,  from  his  earlier  work  on  multi-channel 
time  series  data  analysis  for  use  in  this  study.  Spectral  calculations 
were  done  by  fitting  autoregressive  schemes  to  time  series  data  using 
the  technique  of  a  Fast  Fourier  Transform  (FFT). 


